/VD718981 


AFCRL— 71  — 0012  - 


DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

Cclorado  Stote  University 
Fort  Collins,  Colorado  80521 


MICROWAVE  REENTRY 
PLASMA  DIAGNOSTICS 

8y  James  P.  Rybak 


\ 


Contract  No>  F19628 — 70 — C— 0035 

PROJECT  THEMIS 


Project  No.  7260 


Scientific  Report  No. 3 

CONTRACT  MONITOR:  WALTER  ROTMAN 


January  1971 

MICROWAVE  PKYSICS  LABORATORY 


I  Th'.<  .Hfnrumen}  hos  been  opproved  for  pobiic  releoTe  ond  >^le;  its  distribution  is 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND  UNITED  STATES  AIR  FORCE 
BEDFORD ,  MASSACHUSETTS  01730 


by 

NATIONAL  TECHNICAL 
INFORMj^TlON  SERVICE 

$pr«>sruu!4.  V*  231SI 


AFGRL-7 1-0012 


Depavtiaent  of  Electrical  Engineering 


Colorado  State  University  -  Fort  Collins,  Colorado  80521 


MICROWAVE  REENTRY 
PLASMA  DIAGNOSTICS 

By 

James  P.  Rybak 


v-ontract  No.  F19628-70-C-0035 
Project  No.  7260 


PROJECT  THEMIS 


Scientific  Report  No.  3 
January  1971 


Contract  Monitor?  Walter  Rotman 

Microwave  Physics  Laboratory 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  unllnited. 


Prepared 

for 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COIDIAND 
'JNITED  STATES  AIR  FORCE 


BEDK'RD,  MASSACHUSETTS  01730 


ABSTRACT 

An  analysis  of  surface-mounted  aperture  antenna  admittance  has  been 
performed  in  this  study  to  determine  the  conditions  under  which  admittance 
measurements  can  be  used  for  reentry  plasma  diagnostics.  The  primary 
contribution  of  the  present  work  is  the  determination  chat  the  admittance 
of  a  thin,  microwave  aperture  antenna,  located  on  the  surface  of  a  reentry 
vehicle,  can  be  used  to  obtain  the  values  of  the  electron  density,  electron 
collision  frequency,  ion  sheath  thickness  and  electron  temperature  of  the 
reentry  plasma.  It  is  further  demonstrated,  by  using  admittance  measurements 
made  during  a  reentry  test  f light  as  reported  by  Mayhan  et  al.  (1968  IEEE 
Trans.  Antennas  and  Propagation,  AP-17,  573),  that  open-ended-waveguide 
antenna  admittances  can  be  used  to  determine  the  plasma  electron  density, 
electron  collision  frequency  and  plasma  stand-off  distance  when  the  reentry 
plasma  is  separated,  due  to  aerodynamic  boundary  layer  effects,  from  the 
surface  of  the  reentry  vehicle. 
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CHAPTER  1 


INTRODyCTlON 


The  problems  whicfi  arise  during  the  atmospheric  reentry  of  a  space¬ 
craft  or  missile  are  myriad  and  by  no  means  insignificant.  The  primary 
problem  associated  with  the  reentry  plasma  sheath  is  the  "communii.ations 
blackout"  which  in  many  cases  results  in  a  complete  loss  or  at  least  a 
severe  decrease  in  the  strength  of  radio  frequency  signals  between  the 
reentry  vehicle  and  the  ground.  This  problem  is  of  importance  due  to  the 
loss  of  voice  communicatiu^c  and  data  telemetry  during  the  reentry  of  manned 
space  vehicles  and  the  loss  of  electronic  countermeasures  capability  during 
the  reentry'  of  military  ballistic  missile  payloads  (Jacavanco  1969).  No 
less  important  are  the  concomitant  antee*'^  electrical  breakdown  and  antenna 
pattern  distortion  which  can  also  occuv  before  and  after  the  blackout  period. 
This  problem  lasts  for  up  to  ten  minutes  during  what  is  often  the  most 
crucial  part  of  the  vehicle's  flight.  The  "shuttle"  reentry  vehicies  of  the 
future  will  spend  longer  times  in  the  reentry  phase  of  flight  and  consequently 
will  experience  even  more  prolonged  periods  of  reentry  communications 
problems. 

Considerable  effort  has  been  directed  toward  determining  way.s  to 
alleviate  the  problems  associated  'with  the  reentry  plasma  sheath.  The 
most  promising  techniques  to  lessen  the  reentry  plasma  problem  involve 
the  injection  of  electrophilic  chemicals  into  the  reentry  plasma  flow 
field  upstream  from  the  antenna  locations  and  the  changing  of  the  plasma 
flow  field  by  aerodynamic  shaping  (Huber  and  Sims  196A) .  These  alleviation 
techniques  are  in  the  developmental  stage  and  a  number  of  problems  still 
must  be  solved.  A  thorough  knowledge  of  the  plasma  sheath  properties  is 
needed  in  orde''  to  develop  systems  which  can  maintain  communications  during 
reentry.  Unfortunately,  the  reentry  plasma  sheath  is  the  result  of  chenical 
and  therrocdynamic  processes  'which  are  not  well  undersiond.  In  addition, 
the  .'onditions  which  control  the  rate.s  at  which  these  chemical  and  thermody- 
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nan>ic  processes  occur  change  rapL'Jly  during  reentry,  further  complicating 
the  picture.  Consoquonl ly ,  a  complete  understanding  of  the  plasma  sheath 
cannot  be  obtained  from  -innlytical  approaches.  Sarher,  it  is  necessary 
to  make  accurate  in-tlighl  measurements  of  the  reentry  plasita  sheath 
propertj  es - 

Ideally,  one  would  like  to  deternlne  continuously  the  spatial  variation 
of  the  electron  densitv,  electron  collision  frequency'  and  election  temper¬ 
ature  in  the  inhomogeneous  piasma  sheath  from  Che  surface  of  the  reentry 
vehicle  out  to  the  point  where  the  effect  of  the  plas’na  upon  the  electro¬ 
magnetic  waves  of  interest  becomes  negligible-  However-  the  severe 
environmental  conditions  which  exist  at  ail  but  the  earliest  stages  of 
reentry  preclude  the  use  of  probes  extending  into  the  plasma  to  determine 
the  spatial  variation  of  the  plasma  properties.  Consequently,  it  is  highly 
desirable  to  use  sensors  which  are  flush  mounted  on  the  surface  of  the 
vehicle.  These  sensors  do  not  perturb  the  plasma  flow  field  and  are  not 
subject  to  severe  environmental  conditions.  Unfortunately,  measurements 
from  flush  mounted  sensors  lack  the  spatial  resolution  of  probe  measurements. 
Kevertheless,  surface  mounted  sensor  measurements  are  the  most  practical 
means  of  determining  the  properties  of  the  reentry  plasma  sheath. 

The  radiation  pc  Items  and  admittances  of  surface  mounted  aperture 
antennas  are  changed  markedly  by  the  presence  of  a  reentry  plasma  sheath. 
Consequently,  measurements  of  aperture  antenna  admittance  can  be  used  to 
determine  the  properties  of  the  reentry  plasma  sheath.  The  subject  of 
in-fiight  measurements  v;hich  can  be  related  to  the  properties  of  the 
reentry  piasma  sheath  is  only  in  its  elementary  stages  of  development. 

Much  more  needs  to  be  known  about  the  performance  of  the  surface  mounted 
sensors  in  the  reentry  environment  and  the  proper  interpretation  cf  the 
data  acquired  from  them. 

The  original  work  described  in  this  report  has  been  directed  toward 
furthering  the  knowledge  and  understanding  of  the  behavior  of  surface 
mounted  aperture  antenna>=  in  a  reentry  plasma  environment  and  developing 
techniques  which  enable  one  to  use  aperture  antenna  admittance  measure¬ 
ments  to  determine  directly  the  properties  of  the  reentry  plasma  sheath. 
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These  properties  include  the  electron  density,  electron  collision  frequency, 
electron  temperature  and  the  thickness  of  the  boundary  layer  which  separates 
the  pla&raa  sheath  from  the  reentry  vehicle.  Some  aperture  antenna  con¬ 
figurations  are  able  to  excite  efficiently  longitudinal  or  electroacoustic 
waves  in  addition  to  transverse  waves  in  a  warm  or  compressible  plasma 
fr.e.  a  plasma  with  a  non-zero  electron  temperature).  Consequently, 
the  analysis  of  the  aperture  antenna  sensors  is  based  upon  compressible 
plasma  theory,  where  appropriate,  in  preference  to  the  simpler  but  less 
accurate  incompressible  plasma  theory  which  is  commonly  used. 

The  techniques  developed  in  this  study  are  applied  successfully  to 
measurements  of  aperture  antenna  admittance  made  during  the  reentry  of 
a  Trailbiazer  II  rocket. 

The  equations  appearing  in  this  report  are  written  in  terms  of 
rationalized  mks  units,  exclusively.  Aerospace  and  plasma  physics  jargon 
however  is  not  based  upon  a  consistent  set  of  units.  Data  therefore, 
are  presented  with  commonly  accepted  units. 


CIIAPTKR  li 


APERTUR1-:  ADMITTANCE  FOR  A  RIGID  BOUNDARY 

The  input  admittance  expression  for  a  rectangular  aperture  on  a 
perfectly  conducting  ground  plane  separated  from,  a  compressible  plasma 
half-space  by  an  ion  sheath  of  thickPess  is  fonatilated  as  a  boundar>'- 
value  problem  in  the  current  study.  The  sheath  layer  thickness  is  assumed 
to  be  of  the  order  of  a  Debye  length  (see  Appendix  1} ,  the  sheath  region 
permittivity  is  taken  to  be  that  of  free  space  and  the  sheath  boundary  is 
assumed  to  be  rigid.  The  aperture  admittance  is  evaluated  for  a  number 
of  aperture  dim''nslon3  and  plasma  conditions.  The  plasma  conditions  have 
been  chosen  to  represent,  realistically,  those  encountered  during  reentry. 

These  reentry  plasma  conditions  are  given  in  Table  II-l. 

The  admittance  expression  for  an  open-ended-waveguide  antenna  located 
on  a  ground  plane  covered  with  an  incompressible  (cold)  plasma  has  been 
formulated  and  evaluated  by  Galejs  (1965),  Vllleneuve  (1965),  Compton  (1964), 
Swift  (1967)  and  Mayhan  (1967)-.  A  similar  problem  has  been  treated  by 
Galejs  (1966)  for  a  compressible  plasma  layer.  Galejs  considered  apertures 
id.th  wave-guide  dimensions  only  and  did  not  investigate  the  aperture  admit¬ 
tance  at  frequencies  below  the  plasma  frequency. 

The  determination  of  the  aperture  admittance  expression  in  this  chapter 
is  based  vpon  methods  similar  to  those  used  by  Compton  (1964)  and  Swift  (1967) , 
with  Taodljications  included  to  account  for  the  cooq>re8S'  plasma  effects. 

It  is  fait  that  this  formulation  provides  more  insight  amlng  the  problem 
than  does  the  formulation  developed  by  Galejs  (1966), 

A.  Aperture  Admittance  Expression  The  geixiietry  of  the  problem  is 
shown  in  Figure  IX-1.  The  rectangular  aperture  is  assumed  to  have  a 
prescribed  field  distribution  given  by 

V 

S  «  ~  COB  (ry/b) 

X  & 


Cll-1) 


(II-2) 


for  -a/ 2  <  X  <  all  and  -b/2  <  y  <  b/2  where  V  is  the  aperture 

■  •  ■  o  '  - 

voltage.  The  aperture  field  distributions  given  by  (II-l)  and  (II-2)  are 
valid  for  thin  (a  <  apertures  of  any  length  b  or  for  wid«!r 

ap^ttures  with  b  «  (^lejs  1969).  The  tangential  electric  field 

is  assumed  to  be  identically  zero  on  the  perfectly  conducting  ground  plane. 

^e  effect  of  surface  waves  on  the  aperture  field  distribution  has 
been  neglected  in  this  study.  Galejs  (1966) . h^  shown  that  these 
surface  waves  significantly  affect  the  open-ended-waveguide  antenna  field 
distribution  only  when  the  aperture,  width  is  a  multiple  of  surface  half- 
wavelengths.  Galejs  found  that  this  effect  is  reduced  as  the  plasm 
electron  collision  frequency  is  increased.  The  effect  of  surface  waves  upon 
the  aperture  field  distribution  can  be  assumed  negligible  for  the  values 
of  electron  coli3.sion  . frequency  e^eriehced  during  feehtry. 

The  field  equatd.qns  fbr  the  plasma  regions  are  obtained  by  rewriting 
-  the  linearized  hydrodynamic”  equations  (Oster  1960)  In  the  form 

iujw  H  .1  (II-3) 

o  _  ' 

V  X  H  =  -iwe  E*-  N  ev  (it-4) 

o  o 

(-i«irf-v)mN  V  =  -  N  eE  -  VP  (II-5) 

o  o 

u^mK  V  V  =  iuP 
o 


-K 

V  X  E  = 


(II-6) 
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The  electric  field  E  and  electron  velocity  v  can  be  obtained  in 
terms  of  H  and  F  by  eliminating  either  v  or  E  from  equations  (II-A) 
and  (If-5) .  The  expressions  for  E  and  v  are  then  written 
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It  is  well  known  that  the  vector  components  of  H  can  be  separated 

into  TE  and  Til  parts  which  are  derived  from  two  potential  functions, 

^  and  *  .  The  components  of  H  are  related  to  if*  and  i{»  * 
z  z  z  z 

by 


fs 


1 

H  ^  ^ 

X  Vfi  oy 


H  .  ^ 

y 


2,1.  * 


3 


ibgi^  3x9z 


1 

iuUjj  3y3z 


(Ti-10) 


(II-ll) 


iuu 


f”,s2  I  * 


3y" 


(11-12) 


•V 

Similarly,  the  components  of  E  can  be  written 
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The  potential  functions  ^  snd  T  *  are  those  defined  by 

z  z 

Stratton  (1941)  and  satisfy  the  partial  differentia)  dquationiS 


(11-15) 
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Similarly,  P  satisfies  the  partial  differential  equation  obtained  from 
equations  (ll-o)  and  (11-8) 


+  ~  (uCojrfiv)  -  W  p  «  0  c 


(11-17) 


The  potential  functions  are  defined  both  in  region  I  and  in  region 
II  by  the  relations  (Harrington  1961) 
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r“^(k  ,k  oz)  =  P  ^^(k  ,k  )  exp(ik^^  z) 
X  V  o  X  y  pz 
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The  propagation  coefficients  k  and  k  are  defined  by 
,  sz  pz 
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The  values  of  and  k^^  determined  from  equations  (11-25)  and  (11-26) , 

respectively,  ruust  lie  in  the  first  quadrant  of  the  complex  plane  to  satisfy 

the  radiation  condition  (Spiamerfeld  1952)  at  infinity. 

Equations  (11-18)  through  (IX-24)  are  substituted  in  equation  (TI-8) 

and  sviuations  (11-10)  through  (11-15).  The  boundary  conditions  (Kritz 

Knd  Mintzer  1960)  which  must  be  satisfied  are  1)  continuity  of  E  ,  E  , 

X  y 

H  and  H  at  z  =  0  and  at  z  =  z  ;  and  2)  v  =»  0  at  z  =  z  (rigid 
X  y  o  z  o 

boundary  condition) .  The  mathematics  can  be  simplified  by  defining  the 
following  terms 
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The  tangential  components  of  the  electric  field  must  be  continuous  at 
X  =  0.  Consequently, 


E  =*  -  1 
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k  f(0)F“(2  )  -t-  r— ° - 

y  o'  iioji  c 


o  o 
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(11-38) 
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where  E  and  £  ^  are  the  Fourier  crai^forms  of  the  tangential  electric 
X  y  ' 

field  at  z  =  0. 

liie  tangential,  ^ect^ic  fields  are  identically  zero  on  the  perfectly 
conducting  ground  plane  and  arr  given  by  equations  (II-l)  and  (21-2)  for 
the  aperture  regions.  The  Fourier  trsuisforms  of  the  tangential  electric 
fields  at  2  =  0  are  given  by 
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Eqiiatian  (TIrJ6)  and  the  geor^^try  oi  the  problem  indicate  that 
f.(s)  is  of  t»tc-  form 


f(z)  =  A  exp(ik^  3)  +  B  exp(-ik^  z) 


(11-42) 
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II  II 

The  expressions  for  F  (2  )  and  G  (z  )  now  can  be  detenSined 

o  o 

from  equations  (11-38)  and  (11-39) .  The  results  are 
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with  A,  B,  C,  D  and  defined  by  equations  (11-43),  (11-44),  (11-47), 

(11-48)  and  (11-40),  respectively. 

The  Fourier  transform  of  at  z  =  0  can  be  determined  from  equation 

(II- 11)  along  with  equations  (11-49'  and  (11-50).  The  result  is 


(k  4k  ^)  j  ik  (C-D) 
X  y  >-  ez 


k  (A+B) 
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where  Y  =?  tT"  is  the  characteristic  admittance  of  free  space, 
o  V 

The  time-averaged,  outward  normal,  flow  of  power  per  unit  an 
through  the  aperture  is  given  by  Hessel  et  al.  (1962)  in  the  form 


(11-51) 
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The  pressure  term  P  is  identically  zero  in  the  sheath  region.  Con¬ 
sequently,  equation  (11-52)  can  be  written 


S  = 


(11-53) 
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The  complex  conjugate  of  the  power  radiated  by  the  aperture  can  be 
expressed  using  Parseval's  theorem  (Harrington  1961) 
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The  conjugate  admittance  Y*  of  the  aperture  is  defined  by  the 
relation 
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1  y  2^* 

2  o 


(11-55) 


where  is  the  peak  aperture  voltage. 

Substituting  equation  (11-51)  into  equation  (lT,-54)  and  using 
equation  (11-40)  for  E^^  ,  the  expression  for  the  conjugate  aperture 
admittance  is 
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The  quantity  is  called  the  calculated  aperture  conductance  ana 

3^  is  called  the  calculated  aperture  susceptance. 
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The  evaluation  of  equation  (11-56)  is  simplified  considerably  by 
making  a  transformation  from  rectangular  coordinates  to  polar  coordinates 
by  means  of  the  relations 


S  cos  a  (11-57) 

ky  =  3  sin  a  .  (11-58) 


with  equations  (11-43)  and  (11-44),  defining  A,  B,  C  and  D,  written  in 
terms  of  o  and  0  . 

B.  Evaluation  of  the  Integral  The  integration  of  equation  (11-61) 
cannot  be  done  analytically  and  must  be  done  numerically.  The  method  of 
numerical  integration  chosen  to  evaluate  the  p  integration  in  equation 
(11-59)  is  that  used  by  Richmond  (1966).  This  method  is  said  to  combine 
"the  simplicity  of  the  trapezoidal  rule  v/ith  the  efficiency  of  the  fifth- 
order  Newton-Cotes  formula"  and  has  proven  to  be  quite  satisfactory.  The 
numerical  integration  method  is  outlined  in  Appendix  II.  A  conventional 
Simpson's  rule  procedure  is  used  for  the  a  integration. 

The  integration  of  equation  (11-61)  is  straightforward,  provided  the 
plasma  is  lossy.  For  a  non-lossy  plasma,  surface  wave  poles  and  branch 
points  may  occur  along  the  integration  path.  Seshadri  (1964)  and  Swift 
(1967)  have  described  the  integration  procedure  which  must  be  followed 
when  pules  exist  in  the  path  of  integration.  Only  lossy  plasmas  have 
been  considered  in  the  current  study  to  represent  accurately  the  reentry 
plasma  sheath. 

The  value  of  3  at  which  the  integration  can  be  terminated  depends 

upon  the  dimensions  of  the  aperture.  It  has  been  found  that  the  integration 

can  be  terminated  at  6  =  6.5  for  waveguide  size  apertures  but  that  the 

integration  must  be  carried  out  to  about  B  =  500  fer  very  thin 

(a  <  .Cl  >.  )  apertures.  The  accuracy  of  the  numerical  integration  is 
o 

better  than  one  per  cent. 

C.  Numerical  Calculations  Equation  (11-61)  has  been  evaluated 

numerically  for  a  range  of  aperture  dimensions  and  plasma  parameters. 

Initially,  Che  ion  sheath  thickness  has  been  taken  tc  be  zero.  Figures 

II-2  and  II-3  show  the  aperture  conductance  G  and  susceptance  B  , 

c  c 

normalized  with  respect  to  the  admittance  of  free  space  ,  as  functions 
of  ,  the  real  part  of  the  complex  plasma  dielectric  coefficient.  The 
value  of  the  abcissa  variable  in  these  figures  is  related  to  the  plasma 

frequency  and  to  the  antenna  excitation  frequency  w  through  the 

relation 


normalized  aperture  conductance,  G/y 


Figure  11-3.  Normalized  susceptance  of  aperture  radiating  into 

underdense  plasma;  effect  of  dielectric  coefficient 


It  is  seen  from  figures  II-2  and  II-3  that  plasraa  electron  temperature 
effects,  due  to  the  excitation  of  electroacoustic  waves,  are  small  for 
values  of  greater  than  about  .50  .  Mere  explicitly,  electron  temperature 

effects  are  small  for  values  of  the  antenna  wave  frequency  only  slightly 
greater  than  the  plasma  frequency. 

Also  apparent  from  Figures  I1-2  and  II-3  is  the  fact  that  the  electron 
temperature,  ,  must  be  quite  high  (~10^‘^K)  to  produce  maasureable 
effects  upon  the  aperture  admittance.  The  electron  tempei*ature  at  any 
location  in  the  reentry  plasma  sheath  probably  does  not  exceed  about 
12, 000° K,  This  value  of  electron  temperature  produces  a  negligible  effect 
upon  the  aperture  susceptance  and  only  a  slight  effect  upon  the  aperture 
conductance . 

The  effect  of  the  aperture  length  b  upon  the  calculated  apertv.re 
admittance  is  shov/n  in  Figures  II-A  and  11-5.  The  normalizing  factor 

is  the  free  space  wavelength  of  the  antenna  excitation  signal.  Electron 
temperature  effects  upon  the  aperture  susceptance  become  more  important 
as  the  aperture  length  increases  for  wave  frequencies  above  the  plasma 
frequency.  The  reason  for  this  behavior  is  not  difficult  to  understand 
if  one  examines  equation  (11-61).  The  term 


Cn-63) 

produces  ail  the  electron  temperature  effects  and  is  somewhat  smaller 


than  the  term 


T  =  20000 


ce  cf  apei'tur^  radiating  into 
ffect  of  aperture  length. 
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(11-65) 


becomes  very  small  for  large  values  of  aperture  length  b  except  for 
very  small  values  of  a  ,  This  proc'uces  a  weighting  factor  which  increafes 
the  importance  of  term  (11-63)  with  respect  to  tern  (11-64)  for  larger 
values  of  b  .  Term  (11-64)  contributes  significantly  only  to  the  aperture 
susceptance  while  term  (11-63)  contributes  to  both  the  aperture  conductance 
and  susceptance.  Th-a  result  is  a  linear  relationship  bstwnen  the  apertcro 
conductance  and  the  aperture  length.  .An  examination  of  equation  ((1-61/ 
indicates  that  the  contribution  of  term  (11-63)  to  the  "’•ntegral  is  '•icear 
with  b  due  to  the  a  integration.  Tt’e  electron  temperature  is  observed 
to  determine  the  coefficient  of  proportionality  bet-raen  ths  aperture 
length  and  conductance. 

The  effect  of  the  aperture  length  upon  the  component  of  normalized 
aperture  susceptance  due  to  tho  nonsexo  plasma  electron  temperi;tare 

(tera  IX-63}  is  shown  in  Figure  11-6.  The  value  of  i.as  been  obtained 

by  takiag  the  difference 


8  ’(T  ^0)  =  B  (T  0}  -  B  (T  =0)  .  (11-66) 

c  e  c  K  c  e 

The  value  of  3^’  is  Ouscrved  to  increase  linearly  with  apeiture  length 
b  ,  as  expecteo. 


NORMALlZfO  ELECTROACOUSTIC  SUSCEPTA;^^CE  , 


Figure  II-6.  Nottnalizea  electroacoustic  susct/ptance  of  aperture 

radiating  into  underdense  plasma;  effect  of  aperture 
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ThO  effect  of  the  apetture  width  a  upon  the  normalized  conductance 
and  susceptance  calculated  for  a  sheathless  aperture  excited  at  a  frequency 
slightly  above  the  plasma  frequency  is  shown  in  Figures  II-7  and  II-8, 
respectively.  The  effects  of  the  nonzero  plasma  electron  temperature  are 
diminivhed  as  tha  nov'malized  aperture  width  increases.  Negligible  electron 
temperature  effects  are  observed  when  tha  normalized  aperture  width  is 
.10  or  greater. 

Siectron  tOinparacuve  effects  tend  to  make  the  calculated  aperture 
ausceptance  lej^s  negative  for  very  thin  apertures  but  tend  to  make  the 
aperture  susceptance  more  negative  for  wider  apertures.  This  is  due  to 
the  effect  of  the  weighting  coefficient 

2 

(II-6?) 

upon  terms  (11-63)  and  (T.I-64)  in  equation  (iI-61).  For  large  values  of 
the  aperture  width  a  ,  the  product  of  terms  (11-65)  and  (II-6?)  is  small 
for  ail  but  very  small  values  of  6  .  Plasma  electron  temperature  effects 
in  term  (11-63),  however,  are  most  apparent  for  somewhat  larger  values  of 
8  .  Terms  (11-65)  and  (11-67)  tend  to  suppress  the  electron  temperature 
effects  upon  the  aperture  admittance  as  the  value  of  the  aperture  width 
a  increases. 

The  aC  nittance  of  an  aperture  excited  at  a  frequency  below  the  plasma 
frequency  exhibits  effects  due  to  the  plasma  electron  tempetature  somewhat 
different  from  those  shown  in  Figures  II-2  uhrough  Il-S.  The  calculated, 
normalized  aperture  conductance  and  susceptance  valvies  shc.wii  ir  Figures 
II-9  and  II-IO,  respectively,  are  observed  to  be  msiusureably  affected 
by  plasma  electron  temperatures  encountered  during  reentry. 

Equally  interesting  are  Figures  Il-il  and  11-12  which  show  the  noniaiized 
aperture  conductaftcc  end  susceptance,  respectively,  for  various  values  of 
the  aperture  length  b  ,  The  aperture  is  excited  at  a  frequency  below  the 


NORMALIZED  APERTURE  CONDUCTANCE,  6,/y^ 


NORMALIZED  APERTURE  WIDTH,  o/X. 


Figvire  11-7.  '’onnalized  conductance  aperture  radiating  into 
undeidense  plasmai  effect  of  aperture  width. 
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plasraa  f*.  quency.  The  normalized  aperture  susceptance  and  conductance 
are  observed  Co  vary  linearly  with  the  aperture  length  for  both  zero  and 
nonzero  values  of  the  plasma  electron  temperature. 

The  reasons  for  the  increased  effect  of  the  plasma  electron  temperature 
upon  aperture  admittance  observed  in  Figures  II-9  and  II-IO  together  with 
the  linear  relationship  of  aperture  admittance  to  aperture  IwUgch  observed 
in  Figures  11-11  and  11-12  are  apparent  from  an  examination  of  equation 
(11-61).  The  magnitude  of  term  (11-63)  is  comparable  to  or  larger  than 
the  magnitude  of  term  (11-65)  for  large,  negative  values  of  the  plasma 
permittivity  £  ,  even  for  small  values  of  the  plasma  electron  temperature. 
The  weighting  factor  (I1--64)  along  with  the  sin^  a  and  cos^  a  coefficients 
make  the  contribution  of  term  (11-64)  to  the  integral  in  equation  (11-61) 
negligible  compared  to  Che  contribution  of  term  (11-63),  even  for  relatively 
small  values  of  aperture  length  b  . 

The  effects  of  aperture  width  a  upon  nortaalized  aperture  conductance 
and  susceptance  are  shown  in  Figures  11-13  and  11-14,  respectively.  The 
calculated  aperture  londuccance  is  observed  to  exhibit  smaller  electron 
temperature  effects  a.s  the  aperture  width  is  increased.  Negligible  electron 
temperature  effects  are  observed  for  aperture  widths  greater  than  0.1  X^. 

The  effect  of  Che  plasma  electron  temperature  on  the  normalized, 
calculated  aperture  conductance  and  susceptance  is  shown  in  Figures  II-15 
through  II-17.  It  is  observed  that  temperature  effects  become  more:  pros.viunced 
as  the  ratio  of  the  plasma  frequency  to  the  antenna  excitation  frequ.encv 
becomes  larger.  Temperature  effects,  however,  are  diminished  by  increased 
electron  collisions. 

The  reasons  for  these  effects  can  be  seen  from  an  exarainatio.i  of 
equation  (11-61).  Term  (11-63)  becomes  large  for  smaller  values  of  r> 
as  the  permittivity  of  the  plasma  becomes  more  negative  and  the  electron 
temperature  is  held  constant.  Consequently,  a  given  electron  tempt.iciture 
produces  a  larger  effset  upon  the  admittance  as  the  plasma  permittivity 
becomes  more  negative.  An  increase  in  the  number  of  electron  collisions 
increases  the  imaginary  part  of  the  plasma  permittivity  and  limits  the 
magnitude  whicn  term  (11-63)  can  attain  and,  consequently,  reduces  the 
effect  of  a  given  plasma  electron  temperature  upon  the  apeiiure  edri!  t  L.^ncg . 
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The  previous  admittance  calculationo  hava  been  made  neglecting  sheath 
effects  at  the  surface  of  the  ground  plane.  However,  an  ion  sheath  layer 
will  cover  any  electrically  Isolated  metallic  conductor  in  contact  with  a 
plasma  (Tanenbaum  1967).  The  ion  sheath  thickness  is  several  Debye  lengths 
and  may  appreciably  affect  the  aperture  admittance. 

Figures  II~18  and  11-19  show  the  effect  of  an  increasing  sheath  thick¬ 
ness  upon  the  normalized,  calculated  admittance  of  an  aperture  excited  at 
a  frequency  belov  he  plasma  frequency.  The  iheati.  chickness  is 

normalized  with  respeci  to  the  free  space  antenna  excitation  wave  length. 

The  ion  sheath  has  been  approximated  by  a  free  space  layer  separating 

the  aperture  and  ground  plane  from  the  plasma.  Electron  temperature  effects 

upon  the  admittance  are  seen  to  be  negligible  for  sheath  thicknesses  greater 

than  about  .OCl  X  .  Values  of  the  sheath  thickness  in  excess  of  about 

0 

3  X  10  ^  ,  however,  do  have  an  app.eciabit-  effect  upon  the  aperture 

admittance,  regardless  of  the  plasma  electron  temperature. 

It  must  be  remembered  that,  in  reality,  the  plasma  el  .tron  temperature 
and  the  Debye  sheath  thickness  are  not  independent  quantities.  The  Debye 
distance  X  is  related  to  the  plasma  electron  density  N  and  electron 
temperature  T^  through  the  relation 


o 


(II~68) 


The  admittance  of  an  aperture  excited  at  a  frequency  below  the  plasma 
frequency  and  separated  from  the  compressible  plasma  by  an  ion  sheath  is 
shown  in  Figures  11-20  through  11-23  for  the  range  of  plasma  electron 
temperature  values  normally  encountered  during  reentry.  Both  the  sheath 
thickness  and  the  plasma  electron  temperature  significantly  affect  the 
aperture  admittance  for  the  typical  reentry  conditions  chosen.  The  principal 
effect  of  the  sheath  is  a  large  increase  in  the  aperture  conductance. 
Increased  electron  collisions  are  observed  to  reduce  thi  plasma  electron 
temperature  effects  upon  the  aperture  susceptance.  As  in  the  case  where 
the  sheath  was  neglected,  electron  temperature  effects  are  greater  for 
larger  values  of  („  f  . 

^  p 


Figure  li-18.  Normalized  conductance  of  aperture  radiating  into 
overdense  plasma:  effect  of  ion  sheath  thickness. 
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Noitnalized  susceptanca  of  sheathed  aperture  radiating 
into  overoense  plasma;  effect  of  electron  tenpcrature. 
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Flexure  11-22.  Normalized  conductance  of  sheathed  aperture  radiating  into 
very  overdense  plasma:  effect  of  electron  temperature. 
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Normalized  susceptance  of  sheathed  aperture 
radiating  into  very  overdense  plasma;  effect 
oi  electron  temperature. 


The  conductance  and  susceptance  of  the  sheath  covered  aperture  are 
found  to  vary  linearly  with  the  aportore  length  as  in  the  sheathless  case. 

In  like  'nannar,  plasma  electron  temperature  effects  upor.  the.  admitcan-^e  of 
a  sheath  covered  aperture  are  reduced  as  the  aperture  width  is  increased. 

These  results  are  not  illustrared. 

The  aperture  admittance,  for  a  given  aperture  excitation  frequency, 
is  significantly  at  footed  by  relatively  small  changeo  in  the  plaSina  electron 
density.  This  is  shown  in  figures  11-24  and  11-25  where  calculated  acerture 
conductance  and  susceptance  values,  respectively,  are  shoto  as  functions  of 
frequency  for  tv;o  slightly  different  plasms  layers.  The  curves  labelled 

admittance  values  calculated  for  r  plasma  layer  having  an  electron 
density  equal  to  1.2  x  10 cm  a  collision  frequency  ratio 

v/Wp  equal  to  ,01,  an  electron  temperature  equal  to  12,000'“K 
and  a  sheath  thickness  z  equal  to  three  Debye  le.ngths.  The  curves 
labelled  ^re  admittance  values  calculated  for  a  plasma  layer  having 

an  electron  density  of  2.5  x  10^ ^  cm  ^  with  the  other  parameters  unchanged. 


The  plasma  frequencies  associated  with  the  electron  densities  are  10  and 
14  GHz,  respectively.  It  is  apparent  from  Figures  11-24  and  11-25  that 
significantly  different  aperture  conductances  and  susceptances  occur 
for  a  modest  change  in  plasma  electron  density. 

Figures  11-2  through  11-25  have  sho\m  that,  in  some  cases,  the 
admittance  of  an  aperture  antenna  is  significantly  affected  by  the  electron 


density,  electron  collision  frequency  and  electron  temperature  of  a  con¬ 
tiguous  plasma  layer,  using  the  rigid  boundary  condition.  The  following 
chapter  will  consider  the  aperture  admittance  expression  for  the  case 
when  all  electrons  incident  upor  the  ground  pldue  are  absorbed. 
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CHAPTER  111 

APERTURE  ADMITTANCE  FOR  AIJ  ABSORBING  BOUNDARY 

All  electron^^  coming  within  fp.veral  Debye  lengths  of  uhe  aperturs- 
grouad  plane  surface  were  assumed  repelled  in  the  previous  chapter.  This 
is  not  an  unreasonable  representation  for  a  thin  ion  sheath  and  has  been 
widely  used  in  the  literature.  The  chief  difficulty  with  this  representation 
is  that  the  exact  location  of  the  sheath  edge-  is  not  known  in  iin  actual 
situation  (Watt  1966). 

The  electron  den.'sity  near  a  conducting  boundary  can  be  made  nearly 
uniform  by  eliminating  the  ion  sheath.  The  ion  sheath  can  be  eliminated 
by  making  the  potential  of  the  conducting  boundary  equal  to  the  plaaraa 
space  potential  (Wacserstrom  et  al.  1965).  Kaarlv  ail  the  electrons 
striking  the  conducting  surface  are  then  absorbed.  An  "absorptive  boundary 
condition"  his  been  developed  (Balmain  1965)  using  a  two-sided  velocity 
distribution  in  the  boltzmann  transport  ecuation.  Balmain  used  this  absorptive 
boundaiy  condition  in  his  analysis  of  parallel  plate  and  spherical  RF 
plasma  probes.  The  absorptive  boundary  condition  also  has  been  used  by 
Wait  (1966)  in  ai?  idealized  study  of  the  waves  that  can  propagate  along 
a  p) asma- conductor  interface.  This  bovxndary  condition  generally  has 
not  been  used  in  antenna  analysis. 

A.  Aperture  Admittance  Formula  The  development  of  the  admittance 
formula  proceeds  in  the  same  general  manner  as  before  except  that  nov;  only 
a  single  layer  plasma  geometry  is  considered.  The  potential  functions 
’V  and  T  *  for  this  case  are  written 
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The  propagation  coeffiyient.  k^^  and  k^^  are  defined  in  the  first 
qjjadrant  of  the  complex  plane. 
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The  conpononts  nr  E  and  .1  are  defined  by  equations  (il-iO) 
to  (11-15)  as  before. 

The  boundary  conditions  iaposeu  at  2=0  are  continuity  of 

with  the  absorptive  boundary  condition  (Salraain 

1965;  Wait  1966) 
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inese  ooundary  conditions  be  written  froa  equation  (II-8)  and  equations 
(J.I-10)  through  ai-15)  together  with  equations  (11-62)  through  (11-67)  in 
the  form 


-  (E  -e)  r (k  ‘+k  ^) 

_ JL_ 

N  ec  I  o 
o 


G(o) 


«je  k  P(o) 
+  — 9-Pl- 

K  e 
o 


] 


P(o) 

uK  m 
o 


(III-12) 


E 

X 


I-  kyF(o)  - 


wu  e 
o 


N  e 
o 


(III-13) 


E 

y 


=  0 


k  F(o) 

X 


k  k  G(o) 

_ + 

uy  e 
o 


N  BE 
o 


(111-14) 


The  aperture  field  Fourier  transform  is  defined  by  equation  (11-40). 

The  quantities  F(o),  G(o)  and  P(o)  are  obtained  fivim  equations 
(11-73)  through  (11-75)  in  the  form 
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The  Fourier  transforc  of  at  z  =  0  ca.j  fae  deterained  froo 
equation  (II-ll)  along  with  equations  (11-76)  and  (11-77) .  The  result 
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The  time-averaged,  outward  normal,  flow  of  power  per  unit  area 
through  the  aperture  is  given  by 
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The  conplex  conjugate  of  the  power  radiated  by  the  aperture  can  be 
expressed  using  Parseval’s  theareo  in  the  form 
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The  conjugate  admittance  Y*  is  defined  by 
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Equations  (11-40),  (11-55),  (11-72),  (TI-80)  and  (II-8i)  are  combined 
and  a  change  of  variables  identical  to  that  in  equations  (11-56)  through 
(11-61)  is  performed  with  the  result 
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The  nisnerical  integration  of  equation  {111-21)  is  perfonced  in  the  saine 
sanner  as  '.hat  of  equation  (111-61). 


B,  Hunerical  Results  Aperture  admittance  values  have  been  compr.ted 
for  several  values  of  plasma  parameters  and  aperture  length  for  the  case 
when  the  ion  sheath  is  collapsed  by  a  dc  bias  and  the  absorptive  boundary 
condition  is  valid-  These  admittance  values  are  compared  with  those 
obtained  for  the  case  of  an  ion  sheath  of  zero  thickness  and  the  "hard" 
boundary  condition  in  Figures  lll-i  through  I1I-6. 

It  is  observed  that  the  absorptive  boundary  eliminates  electron 
temperature  effects  upon  aperture  susceptance  and  decreases  electron 
temperature  effects  upon  the  conductance.  Nonzero  electron  temperature 
effects  upon  the  aperture  conductance  are  seen  to  decrease  with  decreasing 
aperture  exciting  frequency.  Figures  III-l  through  III-4  show  that  the 
aperture  admittance  varies  linearly  with  aperture  length  when  the  absorptive 
boundary  condition  is  employed  as  well  as  when  the  hard  boundary  condition 
is  ysed- 

The  aperture  susceptance  is  seen  to  be  constant  over  a  rather  large 
range  of  electron  collision  frequency  values  while  the  zero  temperature 
aperture  conductance  varies  linearly  with  electron  collision  frequency  over 
a  large  range  of  values  when  the  absorptive  boundary  condition  is  employed. 
This  admittance  behavior  is  shown  in  Figures  I1I-5  and  III-6  and  is  similar 
to  that  observed  when  the  hard  boundary  condition  is  employed  for  a  zero 
electron  temperature  plasma. 
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Ths  a<{«Kzt£aace  Cal^J'-catlons  tndicntt!  that  while  eildnatii^jt  thr  1«  '» 
sheath  reduces  the  uncertrtinties  which  exist  concerning  th*?  cxeet  sheath 
thickness,  Che  effect  of  the  nonzero  electron  toaperacare  upon  i.h« 
admittance  is  decreased-  Nevertheless,  aperture  admittance  certsuraaents 
made  at  frequencies  oaly  slightly  below  the  plasma  tre^nency  are  significantly 
affected  by  the  flass^  electron  remperature. 


NORMALJiJeo  APERTURE  CONOU5  TANGE ,  q/y, 


NOP.^JALIZED  APERTURE  LENGTH,  b/'X* 


Figure  Ill-l.  Normalized  conductance  of  aperture  radiating  into 
overdense  plasma;  effect  of  aperture  length  and 
boundary  cenditio/ss. 
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FlRiirc  111-2.  Nor?i,^ll7.ccl  s.uscopca^cc  oi  apfi*-t«trc  caciloting  Into 
overdcn.so  plastna;  effect  of  aperture  length  and 
boundarv  condS  ci&iis. 


WORMALIZEO  aperture  conductance,  G/y 

^  .  tf  o 


Figure  III-3.  Normalized  conductance  of  aperture  radiating  into  very 
ovo.rdense  plasma;  effect  of  aperture  length  and 
boundary  conditions. 
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Figure  III--4.  Normalized  susceptance  of  aperture  radiating  into  very 
overdense  plasma;  effect  of  aperture  length  and  boun¬ 
dary  conditions. 


normalized  collision  frequency,  i/ZtOp 

Normalized  conductance  of  aperture  radiating  into 
overdense  plasma;  effect  of  plasma  losses  and 
boundary  conditions. 
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CHAPTER  IV 


DETERMINATION  OF  REENTRY  PLASJLA  PROPERTIES 
FROM  APERTURE  ADMITTANCE  MEASyREMENTS 

The  calculations  discussed  in  the  previous  chapter  have  shown  the 
effect  of  a  plasma  layer  upon  aperture  antenna  admittance.  The  plasma 
electron  temperatures  experienced  during  reentry  significantly  affect 
aperture  admittance  only  when  the  apei'ture  is  excited  at  a  frequency 
below  the  plasma  frequency  and  only  for  aperture  widths  very  small  com¬ 
pared  to  a  free  space  wavelength-  The  aperture  admittance,  under  these 
conditions,  has  bean  found  to  be  linearly  related  to  the  aperture  length 
for  both  zero  and  nonzero  values  of  the  plasma  electron  ten^jerature . 

These  results  are  in  contrast  to  those  obtained  when  the  aperture  is 
excited  at  a  frequency  above  the  plasma  frequency.  Significant  effects 
of  the  plasma  electron  temperature  upon  the  aperture  admittance  are 
observed  only  when  the  plasma  electron  tSDpcrature  is  considerably 
greater  thar  that  encountered  during  reentry  and  only  when  the  aperture 
length  is  at  least  several  free  space  wavelengths. 

Ion  sheath  effects  upon  the  admittance  of  a  thin  cperture  have 
been  found  to  be  appreciable  and  must  be  included  in  any  realistic 
study.  The  admittance  calculations  also  have  shown  that  the  plasma 
electron  tenperature  has  a  greater  effect  upon  aperture  antenna  admit¬ 
tance  when  electrons  are  reflected  at  the  surface  of  the  ground  plane 
compared  to  that  when  electrons  are  absorbed  by  the  ground  plane. 

Aperture  admittance  vas  calculated  in  Chapters  II  and  III  under 
the  assumption  that  the  plasma  layer  is  infinitely  thick.  In  reality, 
the  thickness  of  a  reentry  plasma  layer  is  at  most  only  of  the  order  of 
centimeters  or  tens  of  centimeters.  However,  plasma  layer  thicknesses 
of  several  centimeters  are  sufficient  to  approximate  an  infinitely  thick 
plasma  when  the  wave  frequency  is  below  the  plasma  frequency.  Consequently, 
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the  antenna-pl"^.aa  layer  ipodei  used  in  the  previous  chapter  is  a  valid 
representation  for  reentry  plasma  layers,  provided  the  antenna  is  excited 
at  a  frequency  below  the  plasma  layer  and  provided  the  plasma  layer  is 
essentially  uniform  over  a  distance  of  several  centimeters. 

Aperture  admittance  measurements  made  for  plasma  layers  not  uniform 
over  distances  of  several  centimeters  require  much  more  sophisticated 
analysis  techniques.  Plasma  layers  of  this  type  often  are  produced  by 
blunted  reentry  vehicles  at  altitudes  above  about  45  km  and  by  sharply 
pointed  reentry  vehicles  above  about  15  km.  For  these  cases,  admittance 
measurements  generally  must  he  made  over  a  wide  raogn  of  frequencies  or  at 
a  number  of  carefully  chosen  discrete  frequencies.  Methods  of  diagnosing 
the  properties  of  grossly  nonuniform  reentry  plasmas  have  been  proposed 
using  admittance  measurement  teamiques  based  on  incompressible  or  cold 
plasma  theory.  It  has  not  been  possible  to  implement  these  methods 
experimentally  (Mayhan  1968b,  1969). 

A*  Thin  Apertures  The  aperture  admittance  expressions  (equations 
(11-61)  and  III-21))  do  not  allow  one  to  determine  the  reentry  plasma 
properties  directly  from  admittance  measurements,  using  a  closed  form 
expression.  However,  it  will  be  shown  in  this  section  that  it  is  possible, 
using  an  iterative  technique,  to  determine  the  plasma  layer  electron 
density,  collision  frequency  and  electron  temperature  from  aperture 
admittance  measurements  made  at  frequencies  below  the  plasma  frequency. 

The  technique  is  not  particularly  sophisticated  and  is  included  merely  to 
show  that  plasma  parameters  can  be  obtained  from  admittance  measurements. 

The  basic  assumption  employed  in  the  following  technique  for  the 
determination  of  plasma  properties  from  aperture  admittance  measurements 
is  that  only  one  value  of  plasma  electron  density  and  electron  collision 
frequency  can  result  in  a  given  value  of  complex  aperture  admittance  if  the 
excitation  frequency  and  electron  temperature  are  specified.  While  this 
assumption  has  not  been  rigorously  proved  true,  it  has  not  been  found  to 
be  uatrue  in  any  cases  investigated  during  the  current  study. 

The  plasma  diagnostic  technique,  based  on  aperture  antenna  admittance 
leeasureraents ,  is  outlined  in  flow  chart  form  in  Figure  IV-1.  Values 
of  aperture  admittance  are  assumed  to  be  known  (measured)  at  two  frequencies 
below  the  plasma  frequen^'y.  Only  one  value  of  admittance  is  considered  at 
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a  time.  An  approximate  value  of  the  plasma  electron  temperature  T^ 

is  chosen,  leaving  the  plasma  electron  density  N  and  collision 

'  o 

frequency  v  required  to  produce  the  measured  admittance  value  with 

specified  at  frequency  u  ,  as  the  unknowns.  The  values  of  and 

\  are  then  calculated  using  the  Newton-Raphson  iterative  technique 
(McCracken  and  Dorn  1964)  to  determine  the  roots  of  the  equations 


G  (N  ,v)  -  G  =  0 

CO  B 


(lV-1) 


(IV-2> 


where  and  3^  are  the  measured  values  of  aperture  conductance  and 

susceptance,  resnectively ,  and  G  and  B  are  the  values  calculated 
■  c  c 

using  equation  (11-61).  The  Newton-Raphson  iterative  technique  is 
discussed  in  Appendix  III. 

Additional  realistic  values  for  T  are  selected  for  the  same 

e 

measured  value  of  admittance  and  the  corresponding  values  of  and 

V  are  calculated.  This  provides  sets  of  values  for  T  ,  N  and 

e  o 

V  ,  all  of  which  corresnond  to  the  measured  value  of  aperture  admittance. 

The  procedure  is  repeated  for  the  admittance  value  measured  at  a 

second  freouencv.  The  result  is  sets  of  T'  N  and  v  values  for  both 

e  o 

measured  admittance  values.  The  set  of  T  ,  H  and  v  values  common 

e  o 

to  both  admittance  values  is  that  of  the  reentry  plasma  layer. 

The  admittance  calculations  performed  in  Chapter  111  have  shown  that 
eliminating  the  ion  sheath  with  a  dc  potential  essentially  eliminates 
electron  temperature  effects  upon  aperture  admittance  in  addition  to 
eliminating  the  uncertainty  associated  with  the  ion  sheath  thickness, 
particularly  when  the  antenna  excitation  frequency  is  only  slightly 
below  the  plasma  frequency.  This  feature  offers  another  possibility 
concerning  the  use  of  aperture  admittance  measurements  for  reentry 
plasma  diagnostics.  This  Second  technique  is  onclined  in  Figure  lV-2. 

The  value  of  aperture  admittance  measured  when  the  ion  sheath  has 
been  eliminated  and  at  a  frequency  about  a  factor  of  ten  below  the  plasma 
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Figure  IV-2.  Flow  chart  of  a  reentry  plasma  diagnostic  method 
based  on  single-frequency  aperture  admittance 


measurements 
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frequency  is  affected  only  slightly  by  the  plesma  electron  temperature. 

This  admittance  measureraeiU  ctn  be  used  to  calculate  the  plasma  layer 
electron  density  and  electron  collision  frequency  v  by  solving 

the  equations 

(IV-3) 

(IV-4) 

using  the  Newton-Raphson  iterative  technique.  The  quantities  and 

are  the  values  of  aperture  conductance  and  susCeptance,  respectively, 
calculated  using  equation  (III-21)  with  u  =  0(T^  =  0).  The  subscript 
m  denotes  measured  values  of  conductance  and  susceptance. 

The  potential  used  to  eliminate  the  ion  sheath  is  then  removed  and 
the  aperture  admittance  is  measured  again,  at  the  same  frequency.  The 
only  unknowns  then  are  the  plasma  electron  temperature  and  ion  sheath 
thickness.  These  quantities  are  then  determined  by  solving  the  set  of 
equations 


(IV-5) 

(XV-6) 


The  values  G  and  5  are  calculated  using  equation  (11-61) . 
c  c 

The  advantage  of  this  technique  is  that  only  single  frequency  adnitt.'<nce 
measurements  need  be  made.  As  far  as  it  is  possible  to  determine  from 
the  open  literature,  the  admittance  of  thin  apertures  never  has  been 
measured  during  reentry.  Consequently,  the  efficacy  of  the  above  mentioned 
techniques  for  determining  the  reentry  plasma  electron  density,  collision 
frequency  and  electron  temperature  cannot  be  demonstrated  using  reentry 


data. 
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The  adcitcance  of  open-ended-waveguit  e  antennas  \.’as  jatasnred  during 
the  reentry  of  a  Trailblazaf  II  rocket  (Caldecott  et  al.  1967).  The 
calculations  in  Chapter  II  have  shown  electron  temperature  effects  upon 
the  admittance  of  apertures  of  this  width  to  be  negiigiale.  Ssowever,  the 
plasma  flow  field  at  the  antenna  locations  was  separated  from  the  aperture 
by  a  distance  ranging  up  to  several  nillinseters  due  to  aerodynanic  boundary 
layer  effects  and  can  be  analyzed  in  a  manner  similar  to  that  outlined 
above. 

The  TraiibJazer  It  flight  on  which  the  opsn-ended-vjaveguide  admit¬ 
tance  was  measured  will  be  described  in  detail  and  the  admittance  date 
analyzed.  First,  however,  the  admittance  characteristics  of  an  open- 
ended-waveguide  antenna  are  discussed . 


B.  Waveguide  Antennas  The  admittance  of  an  open-ended-waveguide 
antenna  can  be  calculated  using  equation  (11-61)  with  appropriate  values 
of  the  aperture  dimensions  a  and  b.  The  value  of  the  electron  acoustic 
velocity  u  is  set  equal  to  zero.  Figures  lV-3  and  IV-4  show  that  plasma 
boundary  layer  thicknesses  or  plarma  ytand-off  distances  of  millimeter 
order  have  an  effect  upon  the  calculated  admittance  of  an  open-ended- 
waveguide  antenna  excited  at  a  frequency  below  the  plasma  frequency.  The 
V7aveguide  antenna  susceptance  is  observed  to  be  virtually  independent  of 
Che  plasma  layer  collision  frequency  for  a  given  value  of  plasma  boundary 
layer  thickness  The  waveguide  conductance,  however,  varies  linearly 

with  Che  plasma  layer  collision  frequency  over  a  large  range.  This 
behavior  of  waveguide  antenna  adraitt..ace  is  analogous  to  that  of  the 
dipole  (King  1961;  Miller  1968a)  and  is  observed  for  frequencies  in  the 

range  1  <  (m  /  m)^  <  6. 

P 

It  is  apparent,  in  this  frequency  range,  that  plasma  layer  electron 
density  and  boundary  layer  thickness  principally  affect  the  waveguide 


antenna  susceptance.  Waveguide  antenna  susceptance  measurements  made  at 


two  frequencies  in  the  range 


1  <  (ai  /u))^ 

P 


determine  th*.  value  of  these  quantities. 


<  6  would  oe  sufficient 
Similarly,  the  waveguide 


to 

antenna 


conductance  is  affected  very  much  by  the  electron  collision  frequcnc>'. 


Measurements  of  the  waveguide  antenna  conductance  can  be  used  to  determine 


the  plasma  electron  collision  frequency,  once  the  other  plasma  parameters 


are  known. 


normalized  waveguide  conductanc 


N0RMAL5ZE0  WAVEGUJDE  SUSCEPTANCE 


The  decenaination  of  the  reentry  plasraa  properties  is  slightly 
more  complicated  if  the  admittance  measurements  are  made  at  frequencies  not 
in  the  range  specified  above.  The  procedure  then  is  similar  to  that 
outlined  initially  (page  63)  foi  thin  apertures.  Figure  T.V-5  shows  a 
block  diagram  representation  of  the  reentry  plasma  diagnostic  technique. 

Values  of  waveguide  antenna  admittance  again  are  assumed  to  be 
known  (measured)  at  two  frequencies  below  the  plasma  frequency.  Ideally, 
the  two  frequencies  should  diff'-r  by  at  least  20  per  cent  of  the  plasma 
frequency.  Only  one  value  of  admittance  is  considered  at  a  t'me.  An 
approximate  valtie  of  the  boundary  layer  thickness  2^  is  chosen.  The 
plasma  electron  density  and  electron  collision  frequency  v  required 

to  produce  the  measured  admittance  value  at  frequency  with  2^  specified 
are  determined  using  the  Newton-Raphson  technique  to  determine  the  roots 
of  the  equations 


G.(N„.v)  ■  ® 


C  o 


m 


B^(N^,V)  -  =  0 


(IV-7) 


(IV-8) 


where  G  ,  G  ,  B  ,  and  B  have  the  same  meanings  as  before.  The  values 
c  m  c  m 

of  G  and  B  are  calculated  using  (11-61)  with  u  =  0(T  =  0). 

c  c  ^ 

Other  realistic  values  of  z  also  are  chosen  for  the  same  admittance 

o 

value  and  the  corresponding  values  of  and  v  are  determined.  This 

provides  sets  of  values  for  z^,  and  v  all  of  which  correspond  to 

the  measured  value  of  waveguide  antenna  admittance. 

The  procedure  is  repeated  for  the  second  admittance  value.  The 

result  is  sets  of  z  ,  N  and  v  values  for  both  measured  admittance 

o  o 

values.  The  set  of  z  ,  N  and  v  values  common  to  both  admittance 

o  o 

values  is  that  of  the  reentry  plasma  layer. 


C.  Trailblazer  II  The  Trailblazer  II  rocket  used  in  reentry- 
test  studies  (Caldecott  et  al.  1967;  Mayhan  et  al.  1968)  is  a  four-stage 
vehicle.  The  first  two  stages  fire  upward  in  the  conventional  manner. 
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Figure  IV-5.  Flow  chart  of  a  reentry  plasma  diagnostic  method 

based  on  two-frequency  waveguide  antenna  admittance 


measurements. 


I 


The  third  and  fourth  stage  engines  face  bactoard  and  drive  the  payload 
package  back  into  the  earth's  atmosphere  along  a  trajectory  nearly 
parallel  to  that  followed  during  the  ascending  part  of  the  flight.  The 
Trsilblazer  II  vehicle  is  capable  of  injecting  a  27  kg  pavload  into  a 
reentry  trajectory  at  a  velocity  of  5.4  kjn  per  second.  The  payload 
reenters  the  atmosphere  about  240  km  frojn  the  launch  point,  minimizing 
the  down-range  support  needed.  The  Trailblazer  II  trajectory  is  shoiTO  in 
Figure  lV-6. 

The  complex  adijiictances  of  two  open-ended  waveguide  antennas  operating 
at  2350  and  5600  Mbiz,  respectively,  v/ere  measured  using  reflectometers 
specially  designed  for  this  purpose  (Bohley  et  al.  1963).  The  antennas 
were  located  at  the  aft  edge  of  the  payload  package  hemispherical  nosecone. 
The  location  of  tae  antennas  with  respect  to  the  stagnation  point  was 
changing  constantly  because  the  payload  was  spin  stabilized  at  13.75 
revolutions  per  second  and  reentered  with  an  initial  attack  angle 
a  of  10.5'*.  Both  antennas  pass.3d  through  positions  A,  B  and  C,  shown 

3 

in  Figure  IV-7,  during  each  revolution,  depending  upon  the  antenna  location. 
The  telemetry  used  to  relay  the  admittance  data  back  to  earth  operated  at 
a  frequency  of  9.21  GHz  and  a  power  of  1100  watts  peak.  The  telemetry 
antennas  were  located  away  from  the  most  severe  plasma  conditions. 

Waveguide  antenna  admittance  was  measured  continually  during  reentry. 

The  reentry'  plasma  was  relatively  thick  and  dense  at  altitudes 
between  about  50  and  30  km.  Consequently,  plasma  properties  az  these 
altitudes  can  be  estimated  from  the  measured  admittance  values,  using 
the  technique  outlined  previously  (page  72).  The  values  of  reentry 
plasma  electron  density  ,  collision  frequency  v  and  boundary  layer 
thickness  z^  ,  determined  in  this  manner  for  several  values  of  reentry 
vehicle  altitude,  are  presented  in  Table  IV-1,  together  with  the  reentry 
vehicle  velocity  and  the  measured  values  of  waveguide  antenna  admittance. 

The  plasma  layer  values  determined  from  the  waveguide  antenna 
admittance  measurements  are  in  agreement  with  those  predicted  by  pre¬ 
flight  air  chemistry  and  fluid  flow  calculations.  The  pre-flight 
calculations  used  siuplified  air  chemsitry  models  and  neglected  the 
reentry  vehicle  angle  of  attacK.  The  predicted  values  of  electron 
density  ,  collision  frequency  v  and  boundary  layer  thickness 

z  were  in  the  ranges 
o 


'A 


ALTITUDE  IN  FEET 


Figure  J.V-6-  Trailblazer  II  flight  trajectory  (Caldecott  et  al- 
1967,  Figure  13). 


Spin  Axis 


figure  VI-1 ,  Hciative  antenna  pcsiticns  during  reentry 
iCaldecoti  ct  a:.  19^7,  Figure  3). 
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and 


5x'0‘^  ^  »viO^-  cm  ’ 

c 

2x10®  <  V  ^  JxiO^  sec"^ 

0  <  2  <  5  tntn 

o 


respectively  (Mayhan  et  al.  l9o8>.  Better  agreement  between  pre-flight 
predictions  and  values  determined  from  admittance  measurements  cannot 
be  expectad  using  simplified  analyses. 

It  is  noteworthy  that  the  pre-flight  predicted  plasma  values  were 
used  by  Mayhan  et  al .  (1968)  to  calculate  the  antenna  admittance  values 
expected  during  reentry.  The  measured  admittance  values  were  in  only 
qualitative  agreement  with  the  predicted  values.  However,  this  agreement 
was  considered  sufficient  to  establish  the  validity  of  the  air  chemistry 
and  fluid  flow  calculations.  No  farther  use  was  made  of  the  measured 
admittance  values  by  Mayhan  et  al. 

The  calculation  of  Trailblazer  II  reentry  plasma  properties  from 
waveguide  antenna  admittance  measurements  has  demonstrated  the  feasibility 
of  using  relatively  simple  techniques  for  reentry  plasma  diagnostics. 
Waveguide  antenna  admittance  measurements,  however,  are  insensitive  to 
the  plasma  electron  temperatures  experienced  during  reentry  and  to  ion 
sheath  or  boundary  layer  thicknesses  less  than  about  a  millimeter^  The 
admittance  of  a  thin  aperture  is  affected  by  these  factors,  making  it 
more  useful  for  reentry  plasma  diagnostics. 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIOxNS  FOR  FUTURE  WORK 

Admittances  of  surface  mounted  aperture  antennas  in  a  reentry 
plasma  environment  have  been  investigated  extensively  in  this  report. 

As  a  result,  a  number  of  significant  conclusions  can  be  stated  and 
several  recommendations  for  future  work  can  be  made. 

Conclusions  The  following  significant  conclusions  can  be 
stated  as  a  result  of  the  current  study; 

1.  The  values  of  e.'ectron  density,  electron  collision  frequency, 
electron  temperature  and/or  stand-off  distance  for  reentry 
plasma  can  be  determined  directly  from  measurements  of  the 
admittance  of  £in  aperture  antenna  located  on  the  surface  of 

a  reentry  vehicle. 

2.  The  excitation  frequency  of  the  aperture  should  be  less  than 
the  plasma  frequency  when  admittance  measurements  are  used  for 
reentry  plasma  diagnostics. 

3.  The  aperture  length  can  be  as  small  as  0.5  of  a  free  space 

RF  wavelength  when  admittance  measurements  are  used  for  reentry 
plasn  diagnostics. 

4.  The  reentry  plasma  electron  temperature  can  be  determined  from 
aperture  antenna  admittance  measurements  if  the  aperture  width 
is  about  0.01,  or  less,  of  a  free  space  RF  wavelength. 

5.  The  effect  of  the  plasma  electron  temperature  upon  aperture 
admittance  depends  upon  whether  the  aperture  ground  plane 
absorbs  or  reflects  incident  electrons, 

6.  The  reentry  plasma  electron  temperature  can  be  determined  from 
aperture  admittance  measurements  if  the  plasma  is  in  contact  with 
the  aperture  ground  plane  or  is  separated  from  the  ground  plane 
by  less  than  0.001  of  a  free  space  RF  wavelength. 
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7.  The  admittance  of  an  open-enaed-waveguide  antenna  can  be  used 
CO  determine  the  plasma  electron  density,  electron  collision 
frequency  and  plasma  stand-off  distance  when  aerodynamic  boundary 
layer  effects  cause  the  reentry  plasma  to  be  separated  from  the 
vehicle  surface.  This  finding  has  been  verified  by  a  companion 
study  at  Colorado  State  University  (Davidson  1970). 

d.  Values  of  waveguide  antenna  admittance  measured  at  two  frequencies 
during  the  reentry  of  a  Traiiblaaer  II  payload  package  were  used, 
in  the  current  study,  to  determine  the  values  of  electron  density, 
electron  collision  frequency  and  plasma  stand-off  dj.stance  asso¬ 
ciated  with  the  reentry  plasma.  These  values  are  cjnsistent 
x-fith  the  values  decerniined  by  analytical  means. 

Recommendations  for  Future  Work  The  analyses  carried  cut  in  the 
current  study  have  been  based  on  the  assumption  of  a  planar  geometry. 

This  is  a  valid  assumption  as  long  as  the  reentry  vehicle  surface  is 
reasonably  flat  over  a  distance  of  several  wavelengths.  The  usefulness 
of  the  current  study  could  be  extended,  however,  by  making  aperture 
admittance  calculations  for  finite  sized  apertures  on  non-planar 
conducting  surfaces.  Diagnostic  measurements  then  could  be  made  at 
locations  on  the  reentry  vehicle,  such  as  the  stagnation  region,  where 
the  surface  definitely  is  non-pianar  and  where  diagnostic  measurements 
could  most  easily  be  related  to  theoretical  predictions. 

The  admittance  calculations  also  have  been  performed  under  the 
assumption  that  the  plasma  layer  is  uniform  and  infinite  in  extent, 
labile  this  assumption  oft  ;n  is  reasonable  when  the  aperture  excitation 
frequency  is  less  than  the  plasma  frequency,  it  would  be  desirable  to 
develop  methods  for  determining  the  spatial  variation  of  the  plasma 
properties  using  aperture  admittance  measurements. 

Mention  has  been  made  of  theoretical  techniques  (Mayban  19*8,  1969). 
based  or-  cold  olasaa  tneory,  vnerein  aperture  admittance  measurements  can 
be  used  to  Jatermine  the  properties  of  a  thin,  inhomogeneous  reentry  plasma 
layer.  These  proposed  techniques  involve  rather  complex  experiiaental  and 
data  analysis  probleais  and  have  not  been  demonstrated  experimentally. 

The  inclusion  of  electron  ceraperature  effects  would  not  make  the  techniques 
significantly  more  complicated  as  long  as  the  plasma  electron  temperature 
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could  be  taken  to  be  approxiraaceiy  constant.  Experimental  admittance 
values  for  this  type  of  nlasma  diagnostics  probably  would  have  to  be 
measured  on  a  swept  frequency  basis  using  techniques  similar  to  those 
described  by  Davidson  et  ai.  (1969). 

The  admittance  values  discus.^ed  previously  all  have  been  calculated 
using  linearized,  hydrodynamic,  compressible  plasma  theory.  This 
hydrodynamic  approach  hsv  been  widely  used  in  the  literature  in  preference 
to  the  kinetic  •'hcory  approach.  Kinetic  plasma  theory  can  provide  a  more 
conplete  and  accvirate  description  of  wave  propagation  in  a  plasma  but  is 
very  difficult  to  apply  to  practical  problems. 

Kinetic  plasma  theory  has  been  used  by  Kuehl  (1966,  1967)  to  calculate 
the  radiation  resistance  of  a  short  filamental  antenna  in  a  warm,  isotropic, 
collisionless  plasma.  He  found  that  the  transverse  component  of  radiation 
resistance  preiicted  fay  kinetic  plasma  theory  agrees  quite  well  with  that 
predicted  by  hydrodynamic  plasma  theory  except  for  frequencies  below  the 
plasma  frequency.  The  longitudinal  component  of  radiation  resistance 
predicted  by  the  two  methods  are  similar  only  for  frequencies  which  are 
large  compared  to  the  plasma  frequency. 

Tlie  analyses  performed  by  Kuehl  involved  many  highly  idealized 
assiimptions  concerning  the  antenna  and  plasoa.  Nevertheless,  Kuehl' s 
results  indicate  the  need  for  more  sophisticated  kinetic  theory  investi¬ 
gations  of  the  performance  of  antennas  in  plasmas.  Antenna  admittance 
values  obtained  using  kinetic  plasma  theory  and  those  obtained  using 
hydrodynamic  plasma  th»  ;ry  should  be  compared  with  experimentally 
3>eosured  values  to  determine  the  conditions  under  which  each  theory 
provides  valid  results. 
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APPENDIX  I 


THE  PLAS'IA  ION  SHEATH 


An  important  concept  in  plasma  physics  is  that  of  the  plasma  ion 
sheath  which  occurs  in  a  plasma  near  a  wall.  The  discussion  of  the  plasma 
ion  shea>:h  in  this  se.ction  is  meant  Co  be  heuristic  and  is  not  mathematically 
rigorous.  A  more  thorough  discussion  of  the  plasma  ion  sheath  is  given 
by  Tancnbaum  (1^67),  The  present  discussion  is  a  summary  of  that  work. 

Charged  particles  incident  upon  a  surface  are,  for  the  most  part, 
lest  to  the  plasma.  The  ions  and  electrons  incident  upon  the  surface 
generally  recombine  and  return  to  the  plasM  as  neutrals.  Also,  some 
electrons  are  able  to  enter  the  conduction  band  of  a  metal  surface 
The  particle  flus  hitting  r%  wall  can  be  calculated  from  a  knowledge  of 
the  velocity  distribution  function  for  the  particles  near  the  wall. 

Assuming  a  one-sided  Maxwellian  velocity  distribution  function  near 
the  wall,  the  electron  flux  and  incident  upon  the  wall  are 
(Tanenbaum  1967) 
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The  electron  flux  to  the  wall  generally  is  one  or  ttro  orders  of 
magnitude  larger  than  the  ion  flux  as  a  result  of  the  large  difference 
between  (T^/m^)  and  (T^/m^)  in  most  plasmas.  Because  of  the  large 
difference  between  the  electron  flu.K  and  the  ion  flux,  an  electrically 
isolated  surface  in  contact  with  a  plasma  develops  a  negati)fc  charge  and, 
hence,  a  negative  potential  with  respect  to  the  plasma.  This  negative  potential 
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provi.des  an  attractive  force  for  the  ions  in  the  plasma  and  a  repulsive 
force  for  the  elections.  The  ion  flux  to  the  wall  increases  and  the 
electron  flux  decreases  until  the  two  fluxes  are  equal  and  the  negative 
potential  at  the  wall  attains  a  fixed  value. 

The  electron  density  in  the  region  near  the  wall  is  lower,  by  about 
one  or  two  orders  of  magnitude,  than  the  electron  density  in  the  main 
body  of  the  plasma.  Furthermore,  near  the  wall  the  electron  density  is 
less  than  the  ion  dens.'ty.  This  region  of  electron  density  and  ion  density 
unbalance  is  referred  to  as  the  plasma  ion  sheath. 

The  thickness  of  the  plasma  ion  sheath  is  of  the  order  of  several 
Debye  distances.  The  Debye  distance  is  defined  by  the  mks  relations 


A  Debye  distance  for  typical  reentry  plasma  parameters  is  of  the 
order  of  10  ®  meter. 
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APPENDIX  IT 


A  NUMERICAL  INTEGRATION  TECHNIQUE 

The  intef;ration  of  equations  (11-61)  and  (111-21)  (pages  16  and 
53,  respectively)  must  be  performed  njmerically,  as  pointed  out  previously. 

A  Simpson’s  rule  technique  is  suitable  for  the  a  integration.  However, 
the  Integrand  varies  ra>  dly  for  some  values  of  B  and  slowly  for  other 
values.  Consequently,  a  numerical  integration  technique  employing  fixed 
intervals  is  not  efficient  for  evaluating  the  6  integration.  Ideally, 
the  intervals  should  be  small  where  the  integrand  varies  rapidly,  and 
should  be  large  ’.Aere  the  integrand  varies  slowly. 

A  computer  program  in  which  the  size  of  the  integration  intervals 
is  adjusted  to  the  \ariation  of  the  integrand  has  been  developed  at 
Ohio  State  University  (Richmond  1965).  The  program  is  said  to 
combine  ’’the  simplicity  of  the  trapezoidal  rule  with  the  efficiency  of 
the  fifth-order  Newton-Cotes  Formula."  The  program  handles  both  real  and 
complex  integrals  of  the  form 

BB 

^  ffx)dx  ’  (All-i) 

AA 

The  computer  program  in  Fortran  IV  is  given  in  Figure  All-1.  The 
variables  ERRR  and  NX  are  the  per  ceni  error  allowed  in  Che  integration 
anti  t*'e  number  of  "gross  intervals"  to  be  used,  respectively.  Only 
the  basic  princinles  of  this  prograia  will  be  discussed. 
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PROGRAM  INT 

COMPLEX  S3S,  SS,  TRAP,  SIl-IP,  FNCP,  TRAZ,  SIMZ,  FNCZ 
READ  (5,1)  AA,  BB,  EPJIR,  NX 
1  FORMA'i  (3E20.3,  13) 

FN  =  NX 

DEL  =  (BB-AA)/FN 
SSS  ^  CMPLX  (0.0, 0.0) 

ERR  =  .01*ERRP./FN 
A  =  AA 

DO  40  NXX  =  1,Na 
IDQC  =  0 
B  =  A  +  DEL 
SS  =  CMPLX  (0.0, 0.0) 

I?X  =  2 
DX  =  DEL/ 2. 

LX  =  1 
X  =  A 
GO  TO  15 
5  TRAZ  =  DX*SS 
MX  -  1 
LX  =  1 
DX  =  DEL 

10  SS  =  CMPLX  (0,0, 0.0) 

LX  =  LX  +  1 
DX  =  .5*DX 
X  =  A  +  DX 
15  DO  20  IX  =  1,MX 
18  SS  =  SS  +  F(X) 

20  X  =  X  +  2.*DX 
IF(LX.EQ.L)  GO  TO  5 
MX  =  2*MX 

21  TRAP  =  .5*TRAZ  +  bX*SS 
DIF  =  CABS(TRAP-TRAZ) 

IF(DIF.GE.r/IP)  MXX  =  IDS  +  1 
DIP  =  DIF 

22  SIMP  =  (4.*TRAP-TRA.Z)/3. 

23  FNCP  =  (16.*SIMP-SIMZ)/15. 

24  ER  =  CABS(1.-FNCZ/FNCP) 

TRAZ  =  TRAP 

SIMZ  =  SIMP 
FNCZ  *-■  FNCP 

25  IF(LX.GT.4)  GO  TO  10 

26  IF(^DCX.GT.4)  GO  TO  30 

27  IF(ER.GT.ERR)  GO  TO  10 
30  SSS  =  SSS  +  FNCP 

40  A  =  A  +  DEL 

“^0  WSITE(6,60)  AA,BB,ERRR,SSS 
60  FORMAT (5E20. 8) 

END 


Figure  All-l. 


A  digital  computer  program  for  numerical 
integration  (Richmond  1966). 
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Let  the  range  of  integration  AA  <  x  ''  BB  be  divided  into  M 
equal  intervsie*  and  let  represent  the  Integrand  value  at  the  end 

points  of  these  intervals  as  in  Figure  AII-2,  Using  the  trapfc=:c'idal 
rule  (HcCracken  and  Dorn  1?64)  -.ith  M  »  1  ,  the  integral  giver,  by 

TRAPj  =>  bX(0.jf.  +  O.Sfj) 

where  DX  =  BB-AA  .  Dividing  the  integration  range  Into  ttwo  ii:terv3ls 
the  trapezoidal  rule  yields 

TKAP^  =  ^  (0.5fj  +  fj  +  O.Sfj) 
while  if  four  intervals  are  used  one  obtains 

TRAPj  “  ^  (Q.5fi  f  fj  +  f j  +  tv  -h  O.Sfj)  . 

Equations  (AI-3)  ana  (AI-4)  can  be  written  in  the  forras 

TRAFj  =  0.5  TRAP,  H-^fj 

nv 

TRAPj  =  0,5  TRAPj  +  -f  (£3  +  f-,)  > 

The  general  relation  is 


TRAP^^j^  =  0.5  T-RAP,  +  (  )  SS 


where  SS  is  the  sunasation  of  the  integrand  values  requii'ed  for 
TRAP.,,  not  included  in  TRAP.  .  Equation  (AlI-7)  appears  i»v  the 

x'tJ.  I 

computer  program  as  statement  2i  where  TRAP  and  TRAZ  denote  TRAP^^^ 


<.\II-2) 


{AII-3) 


(AII-4) 
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and  TRAP^  ,  respectively.  The  quantity  SS  is  evaluated  in  statements 
15  through  20. 

Simpson's  rule  (McCracken  and  Dorn  1964)  yields  the  following  results 
when  used  to  evaluate  the  integral  for  tw’o  and  four  Intervals: 


SL»1P 


DX 


a  =  +  ^^3  +  f5> 


(AII-8) 


~  (4TRAP2  -  TRAPj) 


and 


SIMP,  «  ft  <q  +  +  2£,  +  4£,  +  £,) 


j  (4TRAP3  -  TRAPg) 


(All-9) 


Simpson's  rule  integration  is  seen  to  be  performed  through  a  linear 
combination  of  two  successive  trapezoidal-rule  integrations.  This  is 
carried  out  in  statement  22. 

The  fifth-order  Newton-Cotes  formula  for  numerical  integration 
generally  is  more  efficient  than  the  trapezoidal  rule  or  Simpson's 
rule  and  can  be  obtained  from  a  linear  combination  of  two  successive 
Simpson-rule  integrations  in  the  form 


FNCP^  =  (16  SIMP.  -  SIMP^_j^) 


This  equation  appears  in  the  program  as  statement  23. 


(AlI-10) 
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The  process  of  repeatedly  doubling  the  number  of  intervals  is 
terminated  when  the  results  converge  to  the  specified  accuracy  and  is 
controlled  by  statement  27. 

The  number  of  "gross  intervals,"  denoted  by  NX,  can  have  a  significant 
influence  on  the  time  required  by  the  computer  to  perform  the  integration. 
The  integration  range  (AA  to  SB)  is  divided  into  NX  equal  segments  of 
length  DEL.  The  computer  integrates  over  each  of  these  segments,  one  at 
a  time.  The  value  of  the  total  integral  is  denoted  by  SSS  and  is  obtained 
by  summing  the  subintegrals  (statement  30). 

Intervals  of  equal  length  are  used  over  the  entire  range  from  AA  to 
BB  if  NX  is  unity.  Values  of  NX  larger  than  unity  results  in  the  computer 
using  small  intervals  over  the  segirients  where  the  integrand  varies  rapidly 
and  larger  intervals  over  the  segments  where  the  integiand  varies  less 
rapidly.  The  calculation  time  can  be  minimized  by  a  proper  choice  of 
NX,  although  the  choice  is  not  critical  and  does  not  affect  the  accuracy 
of  the  integration. 
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APPENDIX  ni 

THE  NEWTON-RAPHGON  ITERATIVE  METHOD 

Equations  (IV-1)  and  (IV''2)(page  64)  fonn  a  set  of  two  r.onlinsarj 
ir.te£,rai  equations  in  two  unknowns.  It  is  inpoasible  to  reduos  tha 
problem  sc  that  the  value  of  the  veriabl^ss  and  v  can  be  detencinad 

directly  in  closed  form.  The  Kewton-Raphson  method  (McCracken  and 
Dom  1964)  is  an  iterative  method  for  finding  the  roots  cf  such  equations- 
Let  the  set  of  equations  be  represented  by 


F(x,y>  »  0 


(Alll-i) 


G<x,y)  =  0 


(AIII-2) 


and  let  some  approximate  root.  A  better  approxfm.-'tion  to  the 

root  can  be  obtained  by  writing  a  trin'caCsd  Taylor  series  expansion  of 
the  form 


(AllI-3) 


(AIlI-4) 


The  d-ifiniRg  reistions  for  dx  and  dy  are  obtained  hy  solving 
equations  (AII-3)  and  (AII~4)  and  can  be  written  in  the  form 
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dx  =  - - - - - -  ■ 


vKx  ,y  ) 
n  n 


<AIII-S) 


-92- 


dy  = 


l'\ti  .y  )  -  OU  ,y  ^ 

!i  n  cx  a  r.  n  n  ox  n  n 


(Am-6> 


xijhere 


«  n 


'  V  3G  ^  .  3F  .  3fi 

^  -37  If 


(AIII-7) 


is  assu^^d  rhat  is  not  equal  to  sero- 

A  batter  approximation  to  t'na  yooi  3f  <rquation  (AIII-1)  eni  (AlII-2) 
is  then  given  by 


X  . ,  =  X  +  dx  =  X 
n+1  n  *i 


} 
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T?  ^ 


lx  ,y 
'  n  n 


|-„  3G  „  3F  1 
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+  dy  =  + 


^  3x  ~  ^  3x 
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Successively  better  and  better  approxicsiSions  fce  tuo  root  of  acue«.tiosis 
(AIII-1)  and  (AIIi-2)  are  obtained  by  repeatedly  emp3.o,  Ing  equetions 
(ATII-d)  and  (AIl'1^9).  In  general >  the  sfej-t^cricn  cf  Initial  values 
for  X  and  y  is  not  critical. 
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reentry  plasma.  It  ic  further  demonstrated,  by  using  admittance  measurements 
made  during  a  reentry  test  flight  as  reported  by  Mayhar,  et  al.  (1968  IEEE 
Trsrss.  Antennae  and  i'ropagatlon,  AP-17,  573),  that  open-ended-waveguide 
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plasma  is  esparated,  due  to  aorodynamic  boundary  layer  effects,  from  the 
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